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Abstract 

Heavy quark potentials are investigated at high temperatures. The temperature range 
covered by the analysis extends from T values just below the deconfinement temper- 
ature up to about 4T C in the deconfined phase. We simulated the pure gauge sector 
of QCD on lattices with temporal extents of 4, 6 and 8 with spatial volumes of 32 3 . 
On the smallest lattice a tree level improved action was employed while in the other 
two cases the standard Wilson action was used. Below T c we find a temperature 
dependent logarithmic term contributing to the confinement potential and observe a 
string tension which decreases with rising temperature but retains a finite value at the 
deconfinement transition. Above T c the potential is Debye-screened, however simple 
perturbative predictions do not apply. 
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1 Introduction 



The static quark potential at high temperatures is interesting for several reasons. Phe- 
nomenologically, the properties of quark bound states, in particular of heavy quarko- 
nia, can be derived from potential models quite successfully. It is then important to 
compute the temperature dependence of the potential as this might lead to observable 
consequences in heavy ion collision experiments. Notably, it has been suggested to use 
the suppression of J/^f and production as a signal for the quark-gluon plasma. 
For this purpose, a detailed knowledge of the temperature dependence of the potential 
appears very helpful |2| . 

Moreover, it is well known that a linearly increasing potential at large distances arises 
naturally from a string picture of confinement. As long as one stays in the confined 
phase of QCD, string models then also predict a definite behaviour of the potential at 
finite temperatures p, f|, These predictions ought to be tested by lattice analyses. 

In close vicinity of the deconfinement transition temperature the static quark potential 
and the mass gap i.e. the potential integrated over perpendicular directions are sen- 
sitive to the order of the phase transition. In colour SU(3) the observation of a finite 
mass gap at the critical temperature supported that the transition is of first order 
H, while in SU{2) a continuous decrease to zero with the appropriate Ising critical 
exponents was found @. 

In the deconfined phase, asymptotic freedom suggests that at high temperatures the 
plasma consists of weakly interacting quarks and gluons. Previous numerical studies 
[8-12] have, however, shown that non-perturbative phenomena prevail up to tempera- 
tures of at least several times the critical temperature. In particular, the heavy quark 
potential did not show the simple Debye-screened behaviour anticipated from a re- 
summed lowest-order perturbative treatment fl3| . This might not be too surprising 



as various non-perturbative modes may play a role in the long distance sector of the 
plasma |14] . It is then important to quantify colour screening effects by a genuinely 



non-perturbative approach. 

In the present paper we compute the static quark potential in the pure gluonic sector 
of QCD. We investigate the temperature dependence of the potential over a range of 
temperatures from 0.8 to about 4 times the critical temperature T c . The analysis is 
based on gluon configurations generated on lattices of size 32 3 x N T with iV T = 4, 6 and 
8. This enables us to gain some control over finite lattice spacing artefacts. On the 
smallest lattice a tree level improved gauge action was used while on the two bigger 
lattices a standard Wilson action was employed. We go beyond previous studies of 
the potential in so far the temperature range is covered more densely and also because 



2 



a large set of lattice distances was probed. This helps to extract fit parameters with 
higher reliability. 

The paper is organized such that the next section summarizes theoretical expectations 
on the behaviour of the potential both below and above the transition temperature. 
In section |3| we present and discuss our results for the potential in the confined phase. 
Section |] contains our findings for temperatures above T c and section [| the conclusion. 



2 Theoretical Expectations 



Throughout this paper the potential is computed from Polyakov loop correlations 

(L(0)L\R)) = exp{-V(\R\,T)/T} (1) 



where 



L(x) = -tr Y[ U (x,t) 



(2) 



T = l 



denotes the Polyakov loop at spatial coordinates x. In the limit R — > oo the correlation 
function should approach the cluster value |(L(0))| 2 which vanishes if the potential is 
rising to infinity at large distances (confinement) and which acquires a finite value in 
the deconfined phase. 

In the limit where the flux tube between two static quarks can be considered as a string, 
predictions about the behaviour of the potential are available from computations of 
the leading terms arising in string models. For zero temperature one expects 



V{R) = V -~+*R 



(3) 



where Vq denotes the self energy of the quark lines, a is the string tension and the 
Coulomb-like 1/R term stems from fluctuations of the string |15| . Eq. (§) generally 
gives a good description of the zero temperature ground-state potential although it 
has been shown [|1(J that the excitation spectrum meets string model predictions only 
at large quark pair separations. For non-vanishing temperatures below the critical 
temperature of the transition to deconfinement, a temperature-dependent potential 
has been computed 



as 



V(R,T) 



V 



+ 



7T 1 

— — -arctan(2i?T) 



1 

R 
1 



a T H — T arctanf — — 
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R+^\n(l + (2RTf) (4) 
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In the limit R 3> 1/T this goes over into 



V(R,T) = V + 



IT r 



R + T\n(2RT) (5) 



which had been calculated previously 0]. Note the logarithmic term which originates 
from transverse momentum fluctuations^. So far, it has been left open whether the 
string tension a appearing in eqs. and @ is identical to the zero temperature 
value. In the context of a low temperature or large R expansion, the temperature 
dependent terms appearing in eqs. (§]) and (||) should, however, be considered as 
thermal corrections to the zero temperature string tension. An explicitly temperature- 
dependent string tension was computed by means of a 1/D expansion || 



a{T) 



(7(0) \ I? 



7^2 

1-^ (6) 



where T r was obtained as 



Note, however, that for D — > oo the phase transition is of second order leading to 
a continuous vanishing of the string tension at the deconfinement temperature. In 
colour SU(2), which also exhibits a second order transition, it was established @ that 
cr(T) vanishes ~ (/3 C — (3) u with a critical exponent v taking its 3-D Ising value of 
0.63 as suggested by universality. In the present case of colour SU(3) one expects a 
discontinuous behaviour and a non- vanishing string tension at the critical temperature. 

In the deconfined phase the Polyakov loop acquires a non-zero value. Thus, we can 
normalize the correlation function to the cluster value |(£)| 2 , thereby removing the 
quark-line self energy contributions. Moreover, the quark-antiquark pair can be in 
either a colour singlet or a colour octet state. Since in the plasma phase quarks are 
deconfined the octet contribution does not vanish| and the Polyakov loop correlation 
is a colour-averaged mixture of both 

-V(R,T)/T = }_ -Vi{R,T)/T , -V 8 (R,T)/T (g\ 

9 9 V ; 



At high temperatures, perturbation theory predicts [13| that V\ and V% are related as 

Vx = -8V 8 + 0(g 4 ) (9) 

the context of analyzing numerical data this term has been mentioned in ^| and was 
discussed in detail in 0. 

2 It is, however, small compared to the singlet part. This is true perturbatively, see eq. (|^), as well 
as numerically p0|. 
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Correspondingly, the colour-averaged potential is given by 



V(R,T) _ l lfQB,T) 

T 16 T 2 1 J 

Due to the interaction with the heat bath the gluon acquires a chromo-electric mass 
m e (T) as the IR limit of the vacuum polarization tensor. To lowest order in pertur- 
bation theory, this is obtained as 

M° , (r)V_ /(r) ^ + ^ 



V T J a y ' \ 3 6 

where g(T) denotes the temperature-dependent renormalized coupling, N c is the num- 
ber of colours and Np the number of quark flavours. The electric mass is also known 
in next-to-leading order W71 in which it depends on an anticipated chromo-magnetic 



gluon mass although the magnetic gluon mass itself cannot be calculated perturba- 
tively. Fourier transformation of the gluon propagator leads to the Debye-screened 
Coulomb potential for the singlet channel 

V X {R,T) = -^l e - m ^ R (12) 
R 

where a(T) = g 2 (T)(N^ — l)/(8nN c ) is the renormalized T-dependent fine structure 
constant. It has been stressed |TB[ that eq. ( |T2"D holds only in the IR limit R — > oo 
because momentum dependent contributions to the vacuum polarization tensor have 
been neglected. Moreover, at temperatures just above T c perturbative arguments will 
not apply so that we have chosen to attempt a parametrization of the numerical data 
with the more general ansatz [|l( 



V(R,T) e(T) _ M(T)R 



T (RTY 



e -Kr)H (13) 



with an arbitrary power d of the 1/R term, an arbitrary coefficient e(T) and a simple 
exponential decay determined by a general screening mass fi(T). Only for T ^> T c and 
large distances we expect that d — » 2 and (i(T) — > 2m e (T), eq. (|l0|), corresponding to 
two-gluon exchange. 



3 Results below T c 

The results to be presented here as well as in the next section are based on two different 
sets of data. The first set, referred to as (I) in the following, was generated with a tree- 
level Symanzik-improved gauge action consisting of 1 x 1 and 2x1 loops. The lattice 
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size was 32 3 x 4. We used a pseudo-heatbath algorithm |HJ with FHKP updating [^TJ 
in the SU(2) subgroups. Each heatbath iteration is supplemented by 4 overrelaxation 
steps ||22|| . To improve the signal in calculations of Polyakov loop correlation functions 
link integration ]23| was employed. For each /9-value the data set consists of 20000 to 
30000 measurements separated by one sweep. 

The second set of data (II) was obtained as a by-product of earlier work, the analysis 
of the equation of state [JTSfl . The gauge configurations used in the present study were 
generated with the standard Wilson gauge action on lattices of size 32 3 x 6 and 32 3 x 8. 
The same algorithm as for (I) was employed. The statistics amounts to 1000 to 4000 
measurements separated by 10 sweeps for the N T = 6 data and between 15000 and 
30000 measurements at each sweep for N T = 8. The errors on the potentials as well 
as the fit parameters were determined by jackknife in both cases. 
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Figure 1: The potentials in units set by the temperature at the (3 values analyzed 
for N T = 4 (case I). The critical coupling for this action on our lattice size has been 
determined as (3 C = 4.0729(3) |24|. 

The lattice results for the potential at temperatures below T c are shown in Figures 
[I], ^| and [| The correlation functions, eq. ([!]), have been computed not only for on- 
axis separations but also for some, in the case I almost all, off-axis distance vectors 
R. Although the lattice spacing for the N T = 4 data is larger than for the other two 
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Figure 2: The same as figure |1] except for N T = 6 (case II). The critical coupling is 



(3 C = 5.8938(11) ||. 



lattice sizes rotational symmetry is quite well satisfied due to the use of an improved 
action in this case. As we will focus on the intermediate to large distance behaviour of 
the potential, it was not attempted to specifically treat the deviations from rotational 
invariance at small separations. Note that the distances covered by the data extend 
to RT<4 for (I) while in case II we could obtain signals up to RT < 2. 

The potentials have first been fitted to eq. ([|) with two free parameters, the self- 
energy Vq and a possibly temperature-dependent string tension a{T). These fits work 
rather well even when data at small separations are included because the fit ansatz 
also accounts for a 1/R piece in the potential. The results to be quoted for the string 
tension, Table 1, have, however, been obtained when the data at small separations are 
excluded from the fit. Typically, a minimal distance of RT ~ 1/2 was chosen. The fits 
are stable under variation of -R m i n in this ballpark and return good X 2 values. Varying 
the maximum distance to be fitted does not lead to noticeable changes of the results. 
This holds for all three lattices. 

The results for the string tension, normalized to the critical temperature squared, 
<t/T c 2 , are summarized in Figure |j. The temperature scale has been determined from 
measurements of the string tension at T = |T^[ . The finite temperature string 
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Figure 3: The same as figure |1] except for N T 
P e = 6.0609(9) [M. 



8 (case II). The critical coupling is 



tension is compared to these results at zero temperature, <t(0)/T c 2 , shown as the line in 
the figure. Quite clearly, in the investigated temperature range there are substantial 
deviations from the zero temperature string tension. These deviations amount to 
about 10 % at T/T c = 0.8 and become larger when the temperature is raised. 

Close to T c the results from N T = 6 and 8 at first sight do not seem to agree with the 
numbers coming from the N T = 4 lattices. However, recall that the 577(3) quenched 
theory exhibits a first order transition with the coexistence of hadron and plasma phase 
at the critical temperature. The tunneling rate between the two phases decreases 
exponentially ~ exp(— 2a x (N a /N T ) 2 ) where a = <Jj/Tj} is the normalized interface 
tension. The lattices of data set (II) have a smaller aspect ratio of N a /N T = 4 and 
5.33, respectively, than the 7V T = 4 lattice whose aspect ratio is 8. Correspondingly, 
the ensemble of configurations of the second set contains (more) configurations in the 
"wrong", the deconfined phase. In fact, close to T c , Polyakov loop histograms reveal 
this two-state distribution for N T = 6 and 8, with a clear separability between the two 
Gaussian-like peaks. Such a two-state signal is absent for the N T = 4 data. Carrying 
out the averaging of the potential only over configurations with Polyakov loops in the 
confined peak leads to the corrected data points in Figure 0. At temperatures not 
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Figure 4: The string tension at non- vanishing temperature as obtained from fits with 
eq. (U). The line denotes the zero temperature string tension. In both cases the string 
tension was normalized to T c . For an explanation of the 'corrected' data points see 
text. 



so close to T c this separation of phases is not possible anymore as the Polyakov loop 
histogram has tails into the deconfined phase but it is not clear where one should set 
the cut. 

When we apply this correction the agreement of the results from the three different 
lattices is evident. This shows that the temperature dependence of the string tension 
is not subject to severe discretisation effects. Moreover, the functional form of the 
fit ansatz eq. (f|), as suggested by string model calculations, describes the behaviour 
of the lattice data quite well. However, with increasing temperature we observe a 
substantial decrease of the string tension away from its zero temperature value. Since 
the fit ansatz, eq. contains already a 7r/3 * T 2 term, the decreasing slope of the 
linear part of the potential can not solely be accounted for by this leading correction. 

In order to analyze the linear rising part of the potential in a more model-independent 
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way, in a second round of fits we have compared our data with the ansatz 



V(R,T) = V + a{T)R + CTln{2RT) (14) 

Note that this ansatz differs from eq. ([5]) in so far it summarizes all the linear depen- 
dence on the distance R by an explicitly temperature dependent string tension cr(T). 
Due to its lacking of a 1/R piece this formula is very well capable of describing the 
data but only if the fit is applied to large distances of RT > 1. For data set (II) this 
requirement leaves not too many data points to be fitted. In this case we checked that 



eq. (0) is able to parametrize the potential. However, since we do not have as much 
room to check for stability of the results as one would wish, we refrain from quoting 
results for data set (II). In case I we do have enough distances and obtain fits with 
good X 2 values which are stable under variation of the minimal distance to be included 
in the minimization. 

On data set (I) we clearly observe the logarithmic term contained in eqs. (||) and flljj) . 
The fits return values for the coefficient C of the logarithm which are equal to 1 within 
an error margin of less than 10 %. We thus confirm a logarithmic piece in the potential 
with a strength as anticipated from the string model calculation [|J or, equivalently, 
a subleading power-like l/i?-factor with power 1 contributing to the Polyakov-loop 
correlation function. 

Because of these findings, we fix this coefficient to 1 in the following. The result- 
ing string tension, normalized to its zero temperature value is shown in Figure [5[ 
The temperature dependence compares well with the (modified) prediction || of the 
Nambu-Goto model, eq. (Bj), 



a{T) 



J^2 

1-6^5 (15) 



*(0) T c 

Recall that the string model prediction assumes a second order transition with a con- 
tinuous vanishing of the string tension at the critical temperature. The deconfinement 
transition in pure SU(3) Yang-Mills theory, however, is known to be of first order. 
Thus, a discontinuity at the critical temperature is expected. To account for this, the 
coefficients a and b in eq. fll5|) are allowed to deviate from unity. In fact, the fit to the 
data, shown as the line in figure ||, results in the values a = 1.21(5) and b = 0.990(5). 
This leads to a non-vanishing string tension at the critical temperature of 

^ = 0.121(35) (16) 

This number can be converted into a value for the (physical) mass gap at the transition 
point, m phys (T c )/T c = <j(T c )/T 2 = 0.30(9). This is a bit below but not incompatible 
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Figure 5: The string tension as obtained from fits with eq. fll4|), normalized to its zero 
temperature value. The line is the result of a fit to this ratio with the string model 
motivated ansatz eq. flT5|). The data is compared with the lowest order temperature 
effect on the linear part of the potential, eq. ([|), shown as the dotted line. 



with earlier results of dedicated analyses of the order of the deconfinement transition, 
m-ph ys (T c )/T c = 0.4 — 0.8, as summarized in 0. 

Finally, we compared the string tension cx(T) defined in eq. ( |i~4"D with the leading 
behaviour cr(0) — 7rT 2 /3 as given in eq. (|5|). This is shown as the dotted line in 
Figure |5]. Similarly to Figure |] the comparison fails, reflecting that non-leading terms 
contribute substantially. 
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4 Results above T c 



Above the critical temperature we have normalized the Polyakov loop correlations to 
their cluster value 

(L(O)Lt(i?)) 



V(\R\,T) = -Tin 



(17) 



l(^(0))| 2 

to eliminate the self-energy contributions. In principle, the correlation function itself 
is periodic in R. Alternatively, one can fit the potential, eq. (|I7D , with a periodic 
ansatz, V(R) — > V(R) + V(N a a — R). The second contribution turns out to be very 
small at the distances fitted and both procedures lead to the same results for the fit 
parameters. 
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Figure 6: The potentials V(R, T) for the (3 values above (3 C analyzed on the N T = 4 
lattices. The Polyakov loop correlations have been normalized to their cluster value. 
Potentials and distances are given in units of the temperature. 

In the following we first concentrate on data set (I) which has somewhat better statis- 
tics and which, more important, covers the explored range of distances more densely, 
see Figure |[ As has been explained in section [|, we fit the potentials above T c with 
the generalized screening ansatz, eq. (|TBD, where the exponent d of the Coulomb-like 
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part is treated as a free parameter. It turns out that the value of the exponent and 
the value of the screening mass \i are strongly correlated. In particular at the higher 
temperatures it is difficult to obtain fit results which are stable under the variation 
of the minimum distance included in the fit. These fluctuations have been taken into 
account in our estimates of the error bars. 

At the highest temperatures analyzed we observed that at large quark separations 
the Polyakov-loop correlation decreases below the cluster value. In [18| it was argued 
that finite momentum contributions to the vacuum polarization tensor can give rise 
to a modified screening function which undershoots the exponential Debye decay at 
intermediate distances and approaches the infinite distance limit from below. Despite 
the high, yet limited precision of our data we are not in the position, however, to 
confirm this suggestion. Instead, we have taken an operational approach and have 
added an overall constant to our fit ansatz. 




Figure 7: Fit results for the exponent d of the Coulomb-like part of the potential above 
T c , eq. (|13|), as a function of the temperature. 

In Figure [7] we summarize the results for the exponent. At temperatures very close to 
T c , the exponent d is compatible with 1. When the temperature is increased slightly, 
d starts rising to about 1.4 for temperatures up to 2T C . Between 2 and 3 times T c , the 
exponent centers around 1.5, although the error bars tend to become rather large. A 
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value of 2 as predicted by perturbation theory seems to be ruled out, however, in the 
investigated temperature range. 




Figure 8: Fit results for the screening mass /i(T), eq. (|T3|), as a function of the 
temperature. The lines denote the comparison with the perturbative prediction, 
H(T) = A* m(°)(T), see text. 



The results for the screening mass /x(T) obtained from the same fits with eq. (|T3| ) are 
shown in figure || The screening mass turns out to be small but finite just above T c 
and rises rapidly when the temperature is increased. It reaches a value of about 2.5T 
at temperatures around 1.5T C and seems to stabilize there also. Figure |8] also includes 
a comparison with lowest order perturbation theory, /i(T) = A*m^ > '(T) with m^(T) 
as given in eq. ([TT|). For the temperature dependent renormalized coupling g 2 (T) the 
two-loop formula 



g~ 2 (T)=2b \n 




2 In 



2ttT 



A 



MS , 



was used, where T c / Ajfg = 1.14(4) |26], and the lattice scale was set by the lowest 
Matsubara frequency 2ttT . Perturbation theory predicts the factor A to be 2. Indeed, 
adjusting to the data points at the two highest temperatures, T > 2T C , leads to a 
value of A = 1.82 ± 0.15 which is close to the prediction. However, in view of the 
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results for the exponent d we regard this as an accidental coincidence. This is further 
supported by analyses |T2] in colour SU{2) where the electric gluon mass was obtained 
from gluon propagators and from the singlet potential Vi, see eq (0). Here it was 
found that the observed mass follows a behaviour m e (T) ~ 1.6m^(T). If this result 
could be transferred to the case of SU(3), and [[L0| provides some early evidence for 
this, we would have /x(T) ~ 1 * m e (T) contrary to the perturbative value of 2. 




Figure 9: The screening mass \x from fits with eq. (|T^) with a fixed value of d = 1.5 for 
the exponent. The results from all three different lattice sizes are drawn as a function 
of the temperature. 

The potentials above T c from data set (II) are very similar to the ones already dis- 
cussed. Fits with eq. ( |T3"D with a free exponent do work and return parameter values 
in the same ballpark as in case I. However, because of the much smaller number of 
distances probed in this set the fit results are not as reliable as in case I. Therefore we 
have chosen to carry out fits with eq. (0) but with d kept fixed. For comparison, also 
data set (I) has been treated this way. 

The general feature of these fits is that increasing d from 1.0 to 2.0 leads to decreasing 
numbers for the screening mass. For instance, at 3T C we obtain \xjT ~ 3 for d = 1.0 
whereas with d = 2.0 the result for the screening mass is fi/T m 2. Similar shifts 
occur at all temperatures. The quality of the fits, however, is not always the same. 
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Typically, at temperatures close to T c fits with d = 2 return unacceptable X 2 values 
while for T > 2T C the X 2 values are equally good for all values of d and cannot be 
used to distinguish between the various exponent values anymore. This observation 
fits nicely into the picture as shown in Figure [7]. 

As an example for the temperature dependence of the screening mass at fixed d, in 
Figure |9] we show our results at d — 1.5 for all three different lattice sizes. Recall that 
a value of d ~ 1.5 was favored at all temperatures T>1.2T C of data set (I). The general 
behavior is similar to that shown in Figure §: the screening mass is small close to T c 
and starts to rise quickly. It reaches a kind of plateau with a value of fi/T w 2.5 for 
temperatures between roughly 1.5 and 3 T c . For temperatures beyond 3T C the N T = 8 
data may indicate a slow decrease with rising temperature. The main conclusion to 
be drawn from Figure [5] is that the results from the different lattices, i.e. at different 
lattice spacings are in agreement with each other within the error bars. Thus, in the 
investigated temperature range colour screening effects are not yet properly described 
by simple perturbative predictions. 



5 Conclusion 

In this paper we have analyzed the heavy quark potential at finite temperatures in 
the range 0.8T C up to about 4T C in SU(3) Yang-Mills theory. We have done so on 
lattices with 3 different temporal extents and found results consistent with each other. 
Moreover, the standard Wilson action as well as a tree-level improved Symanzik action 
were used. Again, consistency was observed. This indicates that finite lattice spacing 
artefacts are not futilizing the analysis. 

The potentials at temperatures below the critical temperature of the deconfinement 
transition are well parametrized by formulae which have been derived within string 
models. In particular, the presence of a logarithmic term with the predicted strength 
could be established. However, the obtained string tension shows a substantial temper- 
ature dependence which is not in accord with the leading string model result. Instead, 
we find a decrease of the string tension which is compatible with being proportional 
to (T c — bT) 1 / 2 in the critical region below T c . At the critical temperature the string 
tension retains a finite value of a(T c )/a(0) = 0.121(35), consistent with a first order 
transition. 

Above the deconfinement transition the potentials show a screened power-like beha- 
viour. By comparing the data with perturbative predictions we can further strengthen 
earlier claims that these predictions do not properly describe the potentials up to 
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temperatures of few times the critical one. In particular, it can be excluded that 
the exchange of two gluons with an effective chromo-electric mass is the dominant 
screening mechanism. Judging from the exponent of the 1/R term in the potential, 
at temperatures close to T c it seems that the complex interactions close to the phase 
transition arrange themselves in such a way as to be effectively describable by some 
kind of one-gluon exchange. At temperatures of about 1.5 to 3 times T c we observe 
a behaviour which could be interpreted as a mixture of one- and two-gluon exchange. 
The resulting screening mass scales with the temperature, fi(T) ~ 2.5 T, a pertur- 
bative decrease due to the temperature-dependent renormalized coupling g(T) is not 
really seen. Thus, it is very likely that non-perturbative phenomena and higher order 
perturbative contributions are needed to explain the observed screening behaviour in 
the investigated temperature range. 
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Appendix 



AT R 
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u.y ro 
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u.y ( y 


n 1 11 r >(o\ 
U.lllo^o ) 
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D Q87 


u.uy i u^oj 




n qq^ 
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U.UooZ^lO J 


fi ^ 8n 

O O.oU 


u.ooz 


u.uyi i ^ iu j 


r or 

o.oo 


u.yiy 


U.Uoo4^/o J 


o.oi 


n ore 
u.yoo 


U.UOzO^o J 


5.89 


0.993 


0.0203(13) 






0.0400(20) 


8 6.00 


0.904 


0.0386(6) 


6.04 


0.967 


0.0204(5) 






0.0300(15) 


6.06 


0.9985 


0.0078(7) 






0.0195(15) 



Table 1: Results for the string tension in lattice units from fits with eq. (Q). In those 
cases where two lines are printed for a given (3 the upper line refers to averages over 
the whole sample while the lower line takes into account only configurations which are 
in the confined phase. 
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N T 


T/T c 


aa 2 


4 3.95 


0.804 


0.1592(8) 


4.00 


0.886 


0.1125(4) 


4.02 


0.916 


0.0925(3) 


4.04 


0.947 


0.0686(3) 


4.05 


0.963 


0.0607(3) 


4.06 


0.979 


0.0458(6) 


4.065 


0.987 


0.0316(9) 


4.07 


0.995 


0.0221(11) 



Table 2: Results for the string tension from fits with eq. (]5|). 



20 



